Nonlinear screening in two-dimensional electron gases 
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We have performed self-consistent calculations of the non- 
linear screening of a point charge Z in a two-dimensional elec- 
tron gas using a density functional theory method. We find 
that the screened potential for a, Z — 1 charge supports a 
bound state even in the high density limit where one might ex- 
pect perturbation theory to apply. To explain this behaviour, 
we prove a theorem to show that the results of linear response 
theory are in fact correct even though bound states exist. 

PACS Numbers: 73.20.-r, 73.2Q.Hb 



Screening is a fundamental property of an electron gas 
in arbitrary dimensions. The example of two dimensions 
is of particular interest because of the possible realization 
of quasi- two-dimensional systems in a variety of contexts: 
semiconductor heterostructures ffl , image or band-gap 
surface states at metal surfaces [ff electrons on the sur- 
face of liquid helium , and layered materials . In all 
of these cases, the interaction of external charges with 
the two-dimensional electron gas (2DEG) is a problem of 
both fundamental and practical interest. For example, 
the transport of electrons in a 2DEG is often limited by 
charged impurity scattering and a detailed knowledge of 
the scattering potential is needed for an accurate deter- 
mination of the electron mobility ^ . Scanning tunnelling 
microscopy offers an even more direct means of determin- 
ing the screening response of a quasi-2DEG through the 
observation of adsorbate-induced Friedel oscillations 
Still another class of problems involves the interaction 
with moving charges as might arise in low-energy elec- 
tron scattering Q or tunnelling experiments 0. In this 
case, the dynamic response of the 2DEG is important 
in that electronic excitation, and hence energy loss, will 
occur. 

A charged impurity or projectile typically represents 
a strong perturbation and a nonlinear screening theory 
is in general needed to account for the modifications of 
the local electronic structure. However, in certain situa- 
tions the screened impurity potential may be relatively 
weak and therefore amenable to a perturbative treat- 
ment. This is usually the method adopted to deal with 



donor impurities that are spatially removed from the 
2DEG within a heterostructure j^], although it is rare 
to find quantitative agreement between theory and ex- 
perimentally measured mobilities The situation of 
acceptor impurities within the gas is a much more severe 
perturbation and quite dramatic effects can arise as a re- 
sult of the modified electronic structure |^,|l^. In such 
situations, the screening response has to be determined 
nonlinearly. 

One of our objectives in this Letter is to provide a fully 
self-consistent description of the nonlinear screening in 
an ideal 2DEG within the context of density functional 
theory. A second objective is to use these calculations 
to establish the range of validity of linear response the- 
ory. Somewhat surprisingly, this latter objective is more 
subtle than anticipated as a result of a peculiarity of po- 
tential scattering in 2D, namely the fact that any purely 
attractive potential always has at least one bound state 
pl| , in marked contrast to the situation in 3D. For a pos- 
itively charged impurity we in fact encounter a situation 
in which the screened potential supports a bound state 
even in the high density limit where one would intuitively 
expect a perturbative treatment to be valid. To resolve 
this apparent paradox, we prove what will be referred to 
as the high density screening theorem which states that 
the screening charge can indeed be calculated perturba- 
tively even when bound states exist. In this way we are 
able to justify the use of perturbation theory when at 
first sight it would seem inapplicable. 

The problem we address is the nonlinear screening of 
a stationary point charge, Z, located in the plane of a 
2DEG. The latter is treated as ideal in the sense that 
the electrons are confined to the plane. Of course in real 
applications, such as a heterostructure, a more accurate 
treatment of the electronic states is required. We ignore 
these complications in order to focus on those aspects of 
nonlinear screening which are expected to be independent 
of these details. To stabilize the system, the electrons 
move in the presence of a uniform neutralizing positive 
background. In addition, we assume the electrons to have 
an isotropic effective mass m* and to be immersed in an 
extended dielectric with permittivity e. We use the ef- 
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fective Bohr radius, qq = eh'^/m*e'^, as the unit of length 
and the effective Hartree, H = e^/eao, as the unit of en- 
ergy. The density of the gas, no, is characterized by the 
density parameter = 1/ ^nnQ. 

The static screening response of the 2DEG is deter- 
mined by solving self-consistently the two-dimensional 
Kohn-Sham equations 

- ivV^r) + Ai;cff(r)V'^(r) = E,^/M{r) (1) 



where the effective potential is given by 

Awcff (r) = Ucxt(r) + ^VH{r) + Ai;^c(r) 



(2) 



Here, v^xt = —Z/r is the external potential and Avh 
is the Hartree potential due to the electronic screening 
density, An{r) — n{r) — uq. The change in the exchange- 
correlation potential, Avxci'^c) — Vxcinlj)] — Vxdno], is 
defined in the local density approximation (LDA) using 
the parameterization of the 2D exchange-correlation en- 
ergy given in Ref. . 

The total screening charge is given by 



An{r) ^ J2 I'/^'-WI' + E [I^^WI' - 



(3) 



where the first sum extends over all bound states of the 
effective potential, and the second extends over all oc- 
cupied continuum states up to the Fermi level Ep. We 
assume that each spatial orbital is doubly-occupied for 
spin. The scattering states ipi{r) are shifted asymptot- 
ically relative to the free-particle solutions ipi{r) by a 
phase rjm{E). These scattering phase shifts are related 
to the total screening charge according to the 2D Friedel 
sum rule (FSR) ll| 
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Details of the self-consistent solution of Eqs. (Q) and (||) 
will be presented elsewhere. 

We begin by considering the case of a negatively 
charged impurity {Z = —1), such as an antiproton or ac- 
ceptor state in a semiconductor. Fig. 1 presents the self- 
consistent effective potential Av^s for = 4 (solid curve) 
as a function of the distance from the impurity. This po- 
tential repels electrons almost completely from the impu- 
rity's vicinity, leaving exposed a positively charged disc 
of radius i? ~ which neutralizes the impurity charge. 
As is the case in 3D, this behaviour cannot be reproduced 
in a linear theory. The screened Z — —1 potential does 
not support bound states for any value, as might have 
been expected. We mention this since it has been claimed 
10 that the introduction of a negative test charge into 
a 2D gas can give rise to a potential which is sufficiently 
strong to bind an electron. 

To make contact with this earlier work we compare 
in Fig. |l| the nonlinearly screened potentials with those 



obtained on the basis of linear response theory. The 
chain curve shows the linear response Hartree poten- 
tial (wcxt + Avh) as obtained when local field corrections 
(LFC) are included in the determination of the electron 
screening density |f4|. This potential has a large attrac- 
tive region in real space and supports a bound state for a 
unit negative test charge of one electron mass. This ob- 
servation led to the suggestion of a possible pairing mech- 
anism that could be responsible for a correlation-induced 
instability at low densities [Q. Such a conclusion, how- 
ever, is invalid on two counts. First, the screening of 
the impurity is strongly nonlinear. The dashed curve in 
Fig. |l| shows the corresponding Hartree potential when 
the nonlinear screening density is used. It has a much 
shallower attractive region. But more importantly, an 
electron, as opposed to a negative test charge, also feels 
the effect of the induced xc potential. With this contribu- 
tion included in the full nonlinear potential At^cff (solid 
curve), there is no tendency for bound state formation, 
as confirmed numerically 
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FIG. 1. Comparison of nonlinear and linear screened po- 
tentials: nonlinear with exchange-correlation (solid), nonlin- 
ear Hartree potential (short-dashed), linear Hartree potential 
with LFC (chain); Z = -1, = 4. 

The results for a positive impurity are quite differ- 
ent in that the attractive screened potential supports 
bound states for all densities of physical interest (includ- 
ing ^ 0). For Z = \ there is one m = bound state 
that is doubly occupied. Since the total screening den- 
sity integrates to unity to satisfy the FSR, the continuum 
screening density must itself contribute a total charge of 
+1 in order to compensate for the overscreening provided 
by the bound states. In other words, the Z = 1 impu- 
rity with 2 bound electrons can be viewed as an H^ ion 
which acts as a Z = — I impurity. This was confirmed by 
comparing the Z = \ and Z — —\ continuum screening 
densities. For — 10, the two are virtually the same. 

In Fig. H we show the Z = 1 bound state energy as a 
function of r^. The behaviour seen is surprising in view 
of the corresponding behaviour in 3D. There, the bound 
state energy increases with decreasing since the impu- 
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rity potential is screened more effectively with increasing 
density and, as a result, the bound state eventually ceases 
to exist ll^. Beyond this point, the accuracy of per- 
turbation theory improves with increasing density. The 
contrary behaviour exhibited in Fig. |^ calls into question 
the applicability of perturbation theory in the 2D case. 
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FIG. 2. Bound state eigenvalue vs. r^: nonlinear DFT 
with (solid) and without (dashed) exchange-correlation. 

To address this question, we now prove a theorem re- 
garding electronic screening in the high density limit. 
We consider the introduction of an external potential 
AV^(r) into a uniform noninteracting Fermi gas for ar- 
bitrary dimension D. The parameter A is a coupling 
constant whose physical value is unity. The problem 
is to determine the induced density An(r) due to the 
introduction of XV (r); the effect of interactions will be 
dealt with afterwards. By making use of Dyson's equa- 
tion, G\^s\ = G\ -\- SXGxVGx+sx, for the single particle 
Green function G\{r, r' , z), one can show that the screen- 
ing density satisfies 
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dE 



dY'V{Y') 

X Im [Gx{y,y' ,E + ie)Gx{v\T,E + ie)] 



(5) 



The crucial next step in the argument is to use the an- 
alytic properties of the Green functions to change the 
energy integral from (J^^) to (— /J^)- In the Ep — > oo 
limit it is then permissible to replace Gx by the free- 
particle Green function Go since the energy E can now 
be assumed to be much larger than the strength of the 
(bounded) potential |l^(r)|. Once this is done, the energy 
integration can be changed back to its original range, and 
after integrating with respect to the coupling constant, 
we obtain 

An(r) ~ I dr'V{r') lim ( " dE 

X Im [Go(r, r', E + ie)Go(r', r,E + ^e)] . (6) 

This asymptotic result is valid for any D and applies even 
when the potential V{r) supports bound states. 
An alternative expression for (6) is 



An(r) ^ - hm ^ Xo(q)^(q)e''' ' 
q 



(7) 



which is the result of linear response theory. Here, xo(q) 
is the noninteracting static density response function of 
the system. In the high density limit, q <^ kp for 
all wavevectors for which ^(q) is finite, and we have 
An(r) ~ -xo(0)V(r). In 2D, xo(0) = n'^ and 



An(r) 



-V{r). 



(8) 



Thus we arrive at the interesting conclusion that in 2D 
the screening density takes on a density-independent form 
in the high density limit, and is simply proportional to 
the perturbing potential. 
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FIG. 3. Screening charge density for a model potential: 
bound (dashed) and continuum (chain) state contributions, 
total (sohd). (a) Vo = 0.125 H, (b) V'o = 0.25 H; = 0.5. 

This result can easily be checked numerically. In Fig. ^ 
we show An(r) for a 2D gas with Vs = 0.5 for a model 
potential V{r) — — Vb sin^(27rr/ro)6'(ro — r) which is an 
axially symmetric, double-well potential. With — b 
a.u. and Vb = 0.125 H the potential has a single m — Q 
bound state while for Vb — 0.25 H there are two m — 
bound states. In both cases we see that the total screen- 
ing density is well- approximated by the asymptotic result 
inEq.(|). 

To include the effect of interactions, we identify Awcff 
in Eq. (||) with V{r) and make use of Eq. (||). In the 
high-density limit we then find 
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which in real space gives the Thomas-Fermi potential 
Ip^ . This potential is purely attractive and has a bound 
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state eigenvalue of Eq = —0.2862 H which is the — > 
hmit of the curves in Fig. ^j. This explains why a bound 
state exists in the high density limit and why, in spite 
of this, the result is not in conflict with the applicability 
of linear response theory. It should be emphasized that 
the same argument in 3D leads to the conclusion that no 
bound state can exist in this case. 

As a final practical application we present results of the 
calculation of the energy loss per unit length (or stopping 
power, S) for a projectile of charge Z moving with ve- 
locity V in the plane of the 2DEG. Within the so-called 
kinetic theory framework, the stopping power is given by 
the expression 

S = noVVFcrtriEp) , (10) 

where atriEp) is the momentum-transfer cross section 
defined in terms of the scattering phase shifts by 

4 ^ 

crtriEp) = — ^ sm'^[r]miEF) - rj^+iiEp)] ■ (11) 

m— 

To leading order in the velocity, it is sufficient to deter- 
mine the scattering phase shifts using the static nonlin- 
ear ly screened potentials calculated in the present paper. 



to enhance the stopping power considerably in the range 
— 1 < Z < 1, even to a greater extent than found in 
3D |l^. Finally, comparison with earlier calculations 
|]l9| demonstrates that important differences arise when 
the self-consistently determined nonlinear screening po- 
tentials are used to evaluate the momentum scattering 
cross-sections. 

In summary, we have performed self-consistent calcu- 
lations of the nonlinear screening of a point charge in 
a 2DEG using density functional theory. We have also 
proved a screening theorem which clarifies the behaviour 
of the screening in the high density limit. These results 
find application in a variety of problems, such as charged 
impurity scattering and the stopping power of charged 
projectiles. 
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FIG. 4. Normalized stopping power as a function of 
the projectile charge Z, for Vs — 2: nonlinear screening 
with (solid) and without (dashed) exchange-correlation. The 
straight chain curve gives the quadratic response result. 

In Fig. ^ we show the stopping power as a function of 
the projectile charge Z. For small Z, S has the expansion 
S = SiZ^ + S2Z^ + . . . where the first two terms are the 
linear and quadratic response results, respectively. To 
emphasize the deviations from linear response, we present 
the results in the form S/{vZ'^). In this representation, 
the stopping power including the quadratic response cor- 
rection appears as a straight line with slope S2/V. This 
correction was previously calculated within the quadratic 
random phase approximation and is shown in Fig. ^ 
as the straight line. We can see that corrections beyond 
quadratic response theory are large, especially for pos- 
itive charges. Furthermore, the inclusion of xc is seen 



[1] T. Ando, A. B. Fowler and F. Stern, Rev. Mod. Phys. 

54, 437 (1982). 
[2] P. M. Echenique et al, Chem. Phys. 251, 1 (1999). 
[3] Two-Dimensional Electrons on Cryogenic Substrates, ed. 

by E. Andrei (Kluwer Academic Press, Dordrecht, 1997). 
[4] M. S. Dresselhaus and G. Dresselhaus, Adv. Phys. 51, 1 

(2002); J. Singleton and C. Mielke, Contemp. Phys. 43, 

63 (2002). 

[5] M. F. Crommie et al, Nature 363, 524 (1993); P. A. 
Avouris et al, J. Vac. Sci. Tech. B12, 1447 (1994); P. T. 
Sprunger et al.. Science 275, 1764 (1997). 
[6] T. Nagao et al., Phys. Rev. Lett. 86, 5747 (2001). 
[7] S. Q. Murphy et al, Phys. Rev. B 52, 14825 (1995). 
[8] R. Fletcher et al., Phys. Rev. B 41, 10649 (1990). 
[9] J. Richter et al, Phys. Rev. B 39, 6268 (1989). 
[10] E. Zaremba, Phys. Rev. B 44, 1379 (1991). 
[11] B. Simon, Ann. Phys. (N.Y.) 97, 279 (1976). 
[12] B. Tanatar and D. M. Ceperley, Phys. Rev. B 39, 5005 
(1989). 

[13] F. Stern and W. E. Howard, Phys. Rev. 163, 816 (1967). 
[14] A. GhazaU and A. Gold, Phys. Rev. B 52, 16634 (1995). 
[15] E. Zaremba et al, J. Phys. F 7, 1763 (1977). 
[16] L. Bonig and K. Schonhammer, Phys. Rev. B 39, 7413 
(1989). 

[17] A. Bergara et al, Phys. Rev. B 59, 10145 (1999). 

[18] P.M. Echenique et al, Nucl. lustrum. Methods Phys. 

Res. B 56/57, 345 (1991). 
[19] A. Bret and C. Deutsch, Phys. Rev. E 48, 2994 (1993); 

A. Krakovsky and J. K. Percus, Phys. Rev. B 52, R2305 

(1995); Y.-N. Wang and T.-C. Ma, Phys. Rev. A 55, 2087 

(1997). 



4 



